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1 Introduction 

By the name Euler-Maclaurin, one refers to formulas which relate discrete 
sums to integrals, in particular in the following framework. Let p be a rational 
convex polytope in M. d and let h(x) be a polynomial function on M. d . The 
sum of the values h(x) at integral points of p is written as a sum of terms 
indexed by the set F(p) of faces of p, 

(i) J2 h ^ = E [ D (p>f)-k 

where, for each face f of p , D(p,f) is a differential operator (of infinite order) 
with constant coefficients on M. d . The basic example is the historical Euler- 
Maclaurin summation formula in dimension 1: for a\ < a 2 G Z, 



oi n>l ' n>l 



n! 



«2j 



where 6(n) are the Bernoulli numbers. 

When p is an integral polytope, the existence of such operators is the 
combinatorial counterpart of a homological property of the associated toric 
variety: the invariant cycles corresponding to the faces of p generate the 
equivariant homology module. An Euler-Maclaurin formula amounts to an 
explicit Riemann-Roch theorem, as obtained by Khovanskii and Pukhlikov 
[T7j for integral polytopes corresponding to smooth toric varieties, and ex- 
tended by Cappell and Shaneson [TT] to any integral polytope. Furthermore, 
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by transforming Cappell-Shaneson homological methods into purely combi- 
natorial techniques, valid for any rational (not necessary integral) polytope, 
Brion and Vergne P obtained various expressions for the sum J2 x &pnz d h( x )i 
either as an integral over a deformed polytope followed by differentiation 
with respect to the deformation parameter, or formulas of type (0). 

In this article, we construct differential operators D(p,f), with rational 
coefficients, which satisfy Formula (0) and which moreover enjoy two essential 
properties: they are local and they are computable. The existence of operators 
with these properties was conjectured in [3]. 

By local, one means that D(p, f) depends only on the equivalence class - 
modulo integral translations - of the transverse cone t(p, f) of p along f (see 
Definition |21 and Figure EJ). In particular, if p is an integral polytope, the 
operator D(p, f) depends only on the cone of feasible directions of p along f. 

By computable, one means that there exists an algorithm which computes 
the m lowest order terms of _D(p,f), with running time polynomial with 
respect to the size of the data defining p, at least when the dimension d and 
the number m are fixed. 

On the contrary, Cappell-Shaneson and Brion- Vergne operators are nei- 
ther local nor computable. 

When applied to the constant polynomial h(x) = 1, Formula (0) takes 
the form 

(2) Card(pnZ rf )= £ u (p, f) vol(f) 

fe^(p) 

where the coefficients z/ (p,f) are rational numbers. In the context of toric 
varieties, Danilov [T3j asked the question of existence of coefficients z/ (p,f) 
with the local property, in the case of an integral polytope. This result 
was proven by Morelli jTH] and McMullen [T^]. In Pommersheim and 
Thomas gave a canonical construction of rational coefficients u Q (p,f) which 
satisfy Formula (0), as a consequence of their expression for the Todd class 
of a toric variety. In a companion article, we will similarly obtain a local 
formula for the equi variant Todd class of any toric variety. The result is 
stated in Theorem |2"H 

The computability of our operators D(p, f) extends the following remark- 
able result of Barvinok |2.: when the dimension d is fixed, the number of 
integral points Card(p R Z d ) can be computed by a polynomial time algo- 
rithm. This result is a consequence of Brion's theorem |7], according to which 
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the computation can be distributed over the tangent cones at the vertices, 
and of Barvinok's signed decomposition of a cone into unimodular cones by 
a polynomial time algorithm. Based on this method, efficient software pack- 
ages for integer points counting problems and the effective computation of 
J2xepnz d h( x ) have been developed [T3], (TE], [2HI, [21] • 

Let us now explain the construction of the operators _D(p,f). We define 
the differential operator -D(p,f) through its symbol, using a scalar product 
on M. d . For this purpose, we first associate an analytic function /j,(a) (defined 
on a neighborhood of 0) on R d to any rational affine cone a C M. d . 

For instance, for a half-line s + R+, we have 

e im i 

where [[s]] = n — s, with n6Z,n-l<s<n. Remark that this function 
is analytic at 0, with value at given by fi(s + R + )(0) = ~ — [[s]}. 

The assignment a h- > fi(a) has beautiful geometric properties. The most 
important one is that it is a valuation when the vertex of o is fixed (Theorem 
ED]) . For instance, 

(3) /i(ai U a 2 ) = M a i) + M a 2) - Mai n a 2 ). 

Moreover, /i(a) is unchanged when a is moved by a lattice translation, and 
the map a i— > /x(o) is equivariant with respect to lattice-preserving isometries. 
In particular, with the standard scalar product, we thus get invariants of the 
group 0(n, Z). 

If a contains a straight line, we set //(a) = 0. If a is pointed with vertex 
s, we define /i(a) recursively by the relation: 



(4) Ma)(0 = 




where f denotes a face of a and drrif(x) denotes the canonical Lebesgue mea- 
sure on f defined by the lattice. In (@J), the function /i(t(o, f)) is a priori 
defined only on a subspace of W 1 , namely the orthogonal to the face f. We 
extend it to ~R d by orthogonal projection. Our main point is to show that 
Formula (j3J) actually defines an analytic function. Thanks to the valuation 
property, the proof is reduced to the case of a simplicial unimodular cone. 
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Then we define D(p, f) as the differential operator with symbol /x(t(p, f))(0- 
If the scalar product is rational, the Taylor series of /i(a) has rational coeffi- 
cients, in particular the numbers vo(p, f) in (J2J) are rational. Note that D(p,f) 
involves only differentiation in directions perpendicular to f. 

With this definition, Euler-Maclaurin formula (JIJ for any rational poly- 
tope p follows easily from Brion's theorem. Indeed, the defining formula (@J) 
is formally Formula (JTJ) where the polytope p is replaced by the cone a and 
the polynomial h(x) is replaced by the exponential e^'^ . 

The computability of the functions /i(t(p, f)) is also deduced from Barvi- 
nok's fast decomposition of cones, thanks to the valuation property. 

Moreover, Barvinok proved recently [3] that, given an integer m, there 
exists a polynomial time algorithm which computes the m highest coefficients 
of the Ehrhart quasipolynomial of any rational simplex in IR d , when the di- 
mension d is considered as an input. We hope that our construction of the 
functions fi(a) will lead to another polynomial time algorithm which would 
compute the m highest coefficients of the Ehrhart quasipolynomial for any 
simplex in IR d and any polynomial h(x), when the dimension d and the degree 
of h are considered as input. We want to point out that our construction 
involves only cones of dimension less than m when computing the m highest 
order Ehrhart coefficients. 

In the forthcoming article jT], we will compare our construction to the 
mixed valuation method of j3]. 

Let us illustrate the results in dimension 2. For an affine cone a with 
integral vertex s and edges generated by two integral vectors t>i,t>2 with 
det(t>i,i>2) = 1, (that is to say, a is unimodular), we have: 

/ \ fr\\ 1 , {vi,v 2 ) , 1 1 , 

f(«)(o) = i + — to+(^y)- 

For a general cone, we compute //(a)(0) using the valuation property ©• 

According to Pick's theorem, the number of integral points of an integral 
polygon p C R 2 is given by 

Card(p n Z 2 ) = area(p) + ^length z2 d(p) + 1. 

According to our local formula, as well as Pommersheim-Thomas's, we have: 

Card(p n Z 2 ) = area(p) + -lengthy d(p) + ^ u (p, s) 
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Figure 1: Coefficient z/ (p, s) = /i(t(p, s))(0) at the vertex s 
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Figure 2: Cappell-Shaneson coefficients 



with z/ (p,s) = /i(t(p, s))(0), where s runs over the vertices of p. Thus the 
constant 1 is canonically distributed over the vertices (Figure [T]). Cappell- 
Shaneson coefficients ^2] give a different, non local, distribution of the con- 
stant 1 over the vertices. For instance, the bottom right coefficients in the 
square and the trapezoid of Figure 121 are different, although the vertices have 
the same tangent cone. 

More examples are given at the end of the paper, for polygons p with ra- 
tional (non integral) non unimodular vertices. Based on our Euler-Maclaurin 
operators, we wrote a Maple program which computes the value of the sum 



E, 



epnz 2 x l ^2 



and also the (periodic) coefficients of the corresponding 
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Ehrhart quasipolynomial. 



2 Definitions and notations 

We consider a rational vector space V, that is to say a finite dimensional 
real vector space with a lattice denoted by Ay or simply A. By lattice, we 
mean a discrete additive subgroup of V which generates V as a vector space. 
Hence, a lattice is generated by a basis of the vector space V. A basis of V 
which is a Z-basis of Ay is called an integral basis. 

We will need to consider subspaces and quotient spaces of V, this is why 
we cannot just let V = M. d and A = Z d . The points of A are called integral. 
A point x G V is called rational if qx G A for some integer 9 ^ 0. The 
space of rational points in V is denoted by Vq. A subspace W of V is called 
rational if W n A is a lattice in W. If W is a rational subspace, the image of 
A in V/W is a lattice in V/W, so that V/W is a rational vector space. 

A rational space V, with lattice A, has a canonical Lebesgue measure, for 
which V/A has measure 1. An affine subspace W of V is called rational if 
it is a translate of a rational subspace by a rational element. It is similarly 
provided with a canonical Lebesgue measure. We will sometimes denote 
this measure by drriw- For example, let W be a rational line of the form 
W — s + Mv. Assume that v is a generator of the group M.v fl A (we say that 
v is a primitive vector). Then dmw(s + tv) = dt. 

If Vi G Vq are linearly independent vectors, we denote by □(t>i, . . . ,Vk) 
the semi-open parallelepiped generated by the v^s: 

k 

□(ui, ...,v k ) = y][o, 

1=1 

We denote by vol(D(t>i, . . . , i^)) its relative volume, that is to say its vol- 
ume with respect to the canonical measure on the subspace generated by 
vi, ...,v k . 

We denote by V* the dual space of V. We will denote elements of V 

by latin letters x, y, v , . . . and elements of V* by greek letters £, a, We 

denote the duality bracket by (£, x). 

V* is equipped with the dual lattice A* of A : 

A* = {£ G V* ; (£, x) G Z for all x G A}. 
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If S is a subset of V, we denote by S the subspace of V* orthogonal to 

S: 

S ± = {£ e V* ; (£, x)=0 for all x E S}. 

If W is a subspace of V, the dual space (V/W)* is canonically identified 
with the subspace W L C V*. 

If 5 is a subset of V, we denote by < S > the afline subspace generated 
by S. If S consists of rational points, then < S > is rational. Remark 
that < S > may contain no integral point. We denote by lin(S') the vector 
subspace of V parallel to < S >. 
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Figure 3: The afline space < S > and the linear space lin(S') 

The set of non negative real numbers is denoted by K + . A convex rational 
cone c in V is a closed convex cone Yli=i ~^+ v i which is generated by a finite 
number of elements Vi of Vq. In this article, we simply say cone instead of 
convex rational cone. 

An afnne (rational) cone a is, by definition, the translate of a cone in V 
by an element s e Vq. This cone is uniquely defined by a; it is called the 
cone of directions of a and denoted by dir(o). Thus a = s + dir(a). 

A cone c is called simplicial if it is generated by independent elements of 
Vq. A simplicial cone c is called unimodular if it is generated by independent 
integral vectors v±, . . . , Vk such that vol(D(t> i, . . . , Vk)) = 1. An afline cone a 
is called simplicial (resp. simplicial unimodular) if dir(a) is simplicial (resp. 
simplicial unimodular). 

An afline cone a is called pointed if it does not contain any straight line. 

The set of faces of an afline cone a is denoted by jF(a). If a is pointed, 
then the vertex of a is the unique face of dimension 0, while o is the unique 
face of maximal dimension dim a. 

The dual cone c* of a cone c is the set of £ e V* such that (£, x) > for 
any x G c. 
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A convex rational polyhedron p in V (we will simply say polyhedron) is, 
by definition, the intersection of a finite number of half spaces with boundary 
a rational afline hyperplane. 

Definition 1 We say that p is solid (in V) if < p >= V. 

The set of faces of p is denoted by ^(p) and the set of vertices of p is 
denoted by V(p). 




Figure 4: The transverse cone along an edge in dimension 3 

We now introduce the main geometrical object in our study, the transverse 
cone of a polyhedron p along one of its faces f (see FigureEJ). Let a; be a point 
in the relative interior of f. Recall that the cone of feasible directions of p 
at x is the set c(p, f) := {v G V ; x + eu G p for e > small enough }. It 
does not depend on the choice of x and contains the linear space lin(f). The 
supporting cone of p along f is the afline cone < f > +c(p, f). We denote the 
projection V — > VJ lin(f) by 7Tf. 

Definition 2 Let p be a polyhedron and f a face of p. The transverse cone 
t(p, f) of p along f is the image 7Tf(< f > +c(p,f)) of the supporting cone in 
V/\m(f). 

We will often write simply transverse cone along f, when p is understood. 

If v is a vertex of p, the transverse cone t(p,v) coincides with the sup- 
porting cone v + c(p, v) C V. 
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The transverse cone t(p, f) is a pointed affine cone in the quotient space 
V/lin(f). Its dimension is equal to the codimension of f in < p >. Its vertex 
is the projection 7Tf(x) of any point x of f on V/ lin(f) . 

If a is an affine cone and f is a face of a, then the supporting cone is 
a + lin(f) and the transverse cone t(a, f) along f is just the projection 7r f (a) 
of a on Vyiin(f). 

We shall make use of subdivisions of cones. 

Definition 3 A subdivision of a cone c is a finite collection C of cones in V 
such that: 

(a) The faces of any cone in C are in C. 

(b) If Q 1 and D 2 are two elements of C, then the intersection di fl 2 is a 
face of both Di and $2- 

(c) We have c = U 0G c(). 

If furthermore the elements of C are simplicial cones, the subdivision will 
be called simplicial. 

Example 4 The basic example is the subdivision {R + , (— R + ), {0}} of the 
one- dimensional cone R. 

It is easy to see that any pointed cone admits a subdivision into simplicial 
unimodular cones. 

If S is a subset of V, we denote by x(S) the characteristic function of S 
(also called the indicator function of S). 

Explicit expression of operators D(p, f) involves the Bernoulli polynomials 
b(n,t), defined by the generating series 



The Bernoulli number is b(n, 0) and is denoted by b(n). 

3 Meromorphic functions associated to poly- 
hedra 

By a meromorphic function on V* with rational coefficients, we mean a mero- 
morphic function on the complexification of V* which can be written as the 



(5) 




n=0 
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quotient of two holomorphic functions with rational Taylor coefficients with 
respect to an integral basis of V*. 

We recall the construction of two meromorphic functions with rational 
coefficients on V*, associated to any polyhedron p in V (see the survey [3]). 
The first function J(p) is defined via integration over p, the second function 
S(p) via summation over the set of integral points of p. 

We denote by dm <p> the relative Lebesgue measure on the affine space 
spanned by p. 

Proposition 5 There exists a map I which to every polyhedron p C V asso- 
ciates a meromorphic function with rational coefficients I(p) on V* , so that 
the following properties hold: 

(a) If p contains a straight line, then l(p)=0. 

(b) If £ G V* is such that \e^' x '\ is integrable overp, then 



(d) The map I is a solid valuation: if the characteristic functions x{Pi) 
of a family of polyhedra pi satisfy a linear relation Ylii r iX{pi) = ; then the 
functions I {pi) satisfy the relation 



Example 6 . 

• If p = {s} is a point, then J(p)(£) = e^ ,s K 

• In dimension 1, if p = s + K + , where s E Q, then J(p)(£) = — s j-- 

• If c is a simplicial cone generated by k < d independent vectors V\, . . . , Vk, 
we have 




p 



(c) For every point s G Vq, we have 



I(s + p)(0 = e^I(p)(0. 



r * J (p*) = °- 



{i,<Pi>=V} 



(6) 



I(c)(fl = (-1) 



k vo\(D(v u ...,v k )) 



10 



These formulas follow immediately from the computation in dimension 1, 
J °° e^ x dx = ^ for £ < 0. Thus, for a simplicial cone, the function 7(c) can 
indeed be extended as a rational function on the whole space V*. 

In a similar way, one defines the second meromorphic function, which is 
the discrete analogue of I(p). 

Proposition 7 There exists a map S which to every polyhedron p C V 
associates a meromorphic function with rational coefficients S(p) on V* so 
that the following properties hold: 

(a) If p contains a straight line, then S(p)—0. 

(b) If £ G V* is such that XLepnA < oo ; then 



(7) S(p)(0 = 

rrGpnA 



(c) For every integral point s G A, we have 

(d) The map S is a valuation: if the characteristic functions x(Pi) °f 
a family of polyhedra pi satisfy a linear relation X^ r iX(Pi) — 0, then the 
functions S(pi) satisfy the same relation 

5>S(fc) = 0. 

i 

Example 8 . 

• Ifp = {s} is a point, then we have two cases. If s is an integral point, then 
S({s})(£) = e^ s \ otherwise S({s}) = 0. 

• In dimension 1, if p = s + M + 7 where s e Q, then 



where k is the smallest integer greater or equal than s. 

• Let a be a solid simplicial affine cone with vertex s e Vq. Let vi, . . . , Vd be 
integral generators of the edges of a. Then 



(8) S(a) (0 = 



xe(s+D(v!,...,v d ))nA 



I 11 i _ e (C,f l > 
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We obtain this formula by observing that any element x of the affine cone 
a can be written in a unique way as a sum y + Yli=i n i v i where y lies in 
the semi-open parallelepiped s + . . . , Vd) and the coefficients ni are non 
negative integers, and that the point x is integral if and only if y is. 

Thus indeed S(s + c) can be extended to a meromorphic function on the 
whole of V* . 

Let us check the valuation property in dimension one: 

Example 9 In dimension one, let s 6 Q. The relation x( s + ^+) + x( s ~ 
R + )- X ({s}) =x(K) must imply S(s + R+) + S(s -R+) = S({s}). Ifs = 0, 
we have indeed 

1 +^—, = i = s({s})(0, 



while, ifO<s<l, then 

S(s + R + )(0 



5(a-R+)(0 



1-ef 
1 



l-e-«' 

thus S(s + R + ) + 3(s - R+) = = S({s}). 

The valuation property of the maps p — > S(p) and p — > I(p) have the fol- 
lowing important corollary This was first obtained by Brion [2j using toric 
varieties: 

Theorem 10 (Brion) Let p be a polyhedron in V . Then 

i( P ) = wp> v ))> 

S(p) = s(Mp,v))- 

The singularities of the functions 1(a) and S(a) are easy to compute for a 
pointed affine cone: 

Lemma 11 Let a = s + c be a pointed affine cone with vertex s and let 
Vx, . . . , life be rational generators of the edges of the cone c. The products 

' k \ / k \ 

17 (&«i> )l(s + c)(0 and J] (£,«,) ) S(s + c)(0 
j=i J \i=i J 

are analytic near 0. 
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Proof. It is easy to see that the cone c admits a subdivision into simplicial 
cones whose edges are already edges of c. Thus, thanks to the valuation 
properties of Proposition [7[ it is enough to prove the lemma when c itself is 
a simplicial cone. In this case it follows immediately from Formulas © and 
©. □ 

4 The main construction 

In this section, we will perform the main construction of this article: to any 
affine cone a in V or in a rational quotient V/L of V, we will associate an 
analytic function fi(a) defined in a neighborhood of in V*. In the next 
section, if p is a convex rational polytope in V, we will obtain a local Euler- 
Maclaurin formula for p in terms of the functions /x(t(p, f)) associated to the 
transverse cones t(p, f) of p along its various faces f. 

We will denote the ring of analytic functions with rational coefficients, 
defined in a neighborhood of in V*, by TC(V*) and the ring of meromorphic 
functions with rational coefficients, defined in a neighborhood of in V*, by 
M(V*) 

We will need to extend to the space V* some meromorphic functions which 
are a priori defined only on a subspace of the form (V/L)* = L of V*. For 
that purpose, we fix a scalar product Q(x,y) on V. We assume that Q is 
rational, meaning that Q(x,y) is rational for x,y G Vq. We denote also by 
Q(i,rj) the dual scalar product on V* and we use the orthogonal projection 
proj L ± : V* — > L . If is a meromorphic function (with rational coefficients) 
on L- 1 , we still write <fi for the function on V* defined by £ t— > 0(proj L ±(£)). 
It is meromorphic with rational coefficients. 

Actually, we will do this not only for the space V itself, but also when V 
is replaced by a rational quotient space W. The dual W* is a subspace of 
V*, thus it inherits the scalar product of V*. 

A word of caution. Let L be a rational subspace of V. By means of the 
scalar product Q, we can identify V/L with L ±( 2 C V, the orthogonal of L 
with respect to Q. However these two spaces are not isomorphic as rational 
spaces. The lattice of V/L corresponds to the orthogonal projection of A on 
L ±Q ; it contains the lattice L ±Q fl A, and the inclusion is strict in general, 
see Figure^ 

Let a be an affine cone in V and let f be a face of a. Recall that the 
transverse cone t(o, f) is the projection 7ij(a) of a in V r /Hn(f). When we 
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Figure 5: In dimension 2, the transverse cone along an edge with its lattice 



identify V / lin(f) with the orthogonal lin(f) ±Q , the transverse cone t(a, f) is a 
pointed affine cone in lin(f) ±( 2. 

We denote by C po i n t e d(V) the set of pointed affine cones in V. 

Proposition 12 Let V be a rational space and let Q be a rational scalar 
product on V . There exists a unique family of maps [iw, indexed by the 
rational quotient spaces W of V , such that the family enjoys the following 
properties. 

(a) fi W maps C poin ted{W) to M{W*). 

(b) IfW = {0}, then fi w ({0}) = 1. 

(c) For any pointed affine cone a in W , one has 

(9) S(a) = ^/iin( f) (t(a,f))/(f) 

fe^(a) 

where the sum is over the set of faces of a. 

In Formula @, the transverse cone t(a, f) is a pointed affine cone in the 
quotient space W/lin(f). The function /i W //iin(f)(t( a 5 f)) i s a meromorphic 
function on a neighborhood of in the dual (W/ lin(f))*. We give a meaning 
to the formula by extending this function to a neighborhood of in the whole 
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space W* by means of orthogonal projection. The function 7(f) is defined as 
a meromorphic function on W* as in Section 3. 

Proof. The result is easily obtained by induction on the dimension of 
W. If W — {0}, the only cone is a = {0} and Formula is true. Let a be 
a pointed affine cone in W. Let s be the vertex of a. The transverse cone at 
the zero-dimensional face s is a itself. Formula Q gives 

(10) S(a)(0=e<^W(a)(0+ E WMflW". f)) (0'(f) (0- 

f,dim f>0 

For a face f of positive dimension, the transverse cone t(o, f) is a pointed 
affine cone in the vector space W/lin(f). Therefore, /iwyiin(f) ({•(<*, f)) being 
defined by the induction hypothesis, Formula (fTUj) defines /iw(a) in a unique 
way, as a meromorphic function on 

The following property follows immediately from the definition: 

Proposition 13 If V% C and a is an affine cone contained in V\, then 
the function /iy 2 (a) is the lift to V£ of /ly^a), by the natural restriction map 
V 2 * -> V? . 

In the rest of this article, we will omit the subscript W in the notation /zjy(a). 

Proposition 14 The functions defined in Proposition^^ have the following 
properties: 

(a) For any s E A, one has fi(s + a) = /u(o). 

(b) The map a i— > /x(o) «s equivariant with respect to lattice-preserving 
linear isometries. In other words, let g be a linear isometry of W which 
preserves the lattice A and denote its transpose by t g, then ^(^(a))^ -1 ^) = 

M(a)(0- 

(c) The map a i— > /i(a) multiplicative with respect to orthogonal sums 
of cones. More precisely, if W is an orthogonal sum W = W\ © W2 and 

is an affine cone in Wt for i = 1,2, then 

Mai + 02) = /i(ai)/i(a 2 ). 
^ If a E CpointediW) is such that < a > DA^ = 0, then fiwio) — 0. 

Proof. The invariance in (a) and (b) and the multiplication rule in (c) 
follow immediately from the definition, by induction. To prove (d) , assume 
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that < a > flAy = 0. Then 5(a) = and, for any face f of a, the transverse 
cone t(o, f) does not contain any integral point of Vyiin(f), therefore, by 
induction, /i(t(a, f)) = for dimf > 0, hence //(a) = 0. □ 

If a is an afhne cone in V which contains a straight line, we define /x(o) = 
0. Since all faces f of a contain a straight line, Formula © still holds in this 
case. Thus we have defined //(a) for any rational affine cone in any quotient 
space of V. 

Our objective is to show that fi(a) is indeed analytic near 0, but we will 
first observe some further properties of this family of functions. 

It is easy to compute // in dimension 1. Let t G Q such that < t < 1. 
Then //(-* + K+)(£) is defined by 

1 f°° 

S(-t + R+) (0 = 7 = e~VH + K +) (0 + / ^dx 

1 — e J-t 

hence 

(ii) mH + i + )(0 = t ^ + |. 

We may write this in terms of the Bernoulli polynomials b(n,t), defined by 
the generating series (JHJ). We obtain 



d2) (i( _ t+R+ ) K ) = _ E K(^±M r . 

n=0 



71+ 1) 



Let a be a 1-dimensional pointed cone in V, that is to say a half-line. If 
a does not contain any integral point, then //(a) =0. If a contains integral 
points, then there exists an integral point a e a such that the translated 
half-line a — a is of the form 

(-t + R+)v 

where v is a primitive integral vector and t G [0, 1[. 
By a similar computation, we get, for £ G V*, 

(13) Ma)(0 = ^-77^ + 



The next step is crucial to our construction; we will prove that the map 
c — > fi(s + c) enjoys the valuation property: 
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Proposition 15 Let Cj be a finite family of cones in V . Assume that there 
exists a linear relation between their characteristic functions 5^i r «x( c i) = 0- 
Then, for any s 6 Vq, we have the corresponding relation £V r iA t ( s + c «) = 0. 

Example 16 Consider the subdivision of Example^ in dimension one. Let 
s E Q. Then 

= = fx(s + R + ) + //(s - R + ) - //({s}). 
Indeed, if s is an integer, we have 

//(s+R+XO+Ms-R+XC) = (73^ + ^) = 1 = 

while, if s e] — 1, 0[, we have 

A*(«+R+)(0+A*(a-K+)(0 = C^+j) +(rr^-^) = = MW)(0- 

Proof. 

We will prove the proposition by induction on dim\A By a standard 
argument ([2T]. see also [S]), it suffices to prove the result in the following 
particular case (see Figure EJ). Let a be a pointed solid affine cone in V with 
vertex s, let H be an affine hyperplane through s. Denote by H ± the closed 
half-spaces separated by H. Then we have 

X (a) = x(o n #+) + x(a n #~) - x(a n #) 

and we must prove: 

(14) n(a) = /i(o n H + ) + //(a n H~) - fi(a n iZ). 
We proceed to prove (JHJ). 

The functions 5"(a) have the valuation property 

5(a) - S{a n #+) - 5(a n H~) + 5(o n H) = 0. 

By applying Formula we obtain the following expansion of the left hand 
side 

(15) £>(l(a,f))I(f)- £ Mt(anif + ,f))/(f) 

fe^(a) fe^(anH+) 

- ]T //(t(onF-,f))/(f)+ ]T /i(t(anF,f))/(f). 
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Figure 6: The cut of a 3-dimensional cone with 5 edges by a 2-dimensional 
plane H. The cone is represented by its slice in the figure plane. Thus H is 
represented by a line. 



Let L be a affine subspace of V of dimension > 0. We will show 

(16) Yl Mt(a,f))/(f)- £ Mt(aniJ+,f))/(f) 

feJ r (a),<f>=L fej r (an/i'+),<f>=L 

//(t(anir,f))i(f)+ ]T Mt(an^,f))/(f) 

fe.F(anH-),<f>=z fe.F(an.ff),<f>=£ 

= 0. 

From the relation (|16|h it follows that the terms in (|15|l corresponding 
to the faces f of positive dimension add up to 0. Therefore the contribution 
of the 0-dimensional vertex {s} to (|T5|) is also equal to 0, which proves the 
relation (JUJ). 

We fix L and we proceed to prove 

Remark that all the transverse cones which appear in (|T6*j) are affine cones 
in vyiin(L). We will apply the induction hypothesis to V/\m(L). 

I) First we consider the case where there is a face f of a such that < f >= 
L. There are three cases, according to whether the relative interior of f meets 
both the interiors of H ± , or only one, or none of them (in the third case, f is 
contained in H). 

• Case 1.1: The relative interior of f meets both the interiors of H ± . Then 
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Figure 7: Case 1.1 



f fl if ± is a face of a PI if ± and < f n >= L. Thus we have to prove 

Mt(a, f))/(f)-Mt(an# + , fnif + ))f (fnif + )-//(t(anff-, fDH-))I(fr)H-) = 0. 

The three transverse cones t(o, f), t(a n if + , f n if + ) and t(o n if - , f n if ~) 
coincide. The integrals add up: 

f(f) = f(fnff+) + f(fnff-), 

thus we get 

/i(t(a,f))(/(f)-/(fnii + )-/(fnii-)) 

which is equal to as required. 

• Case 1.2: The relative interior of f is contained in the interior of, say, 
if + . Then f = L fl a fl H + is also a face of a n H + , but L fl a fl H and 
L fl a fl H are smaller dimensional, or empty. This time we have to prove 

Mt(a, f))i(f) - /x(t(a n ii + , f n ii+))/(f n /f + ) = o. 

The transverse cones t(o, f) and t(o fl H + , f) coincide, so we get 

Mt(a,f))(/(f)-/(fnii + )). 

As f = f fl H + , this is equal to as required. 

• Case 1.3: f is contained in if, thus it is a face of all four cones, in other 
words 

f = f n H + = f n H~ = f n H. 



19 



hspacelcm 
Figure 8: Case 1.2 

/ <f>=L 




Figure 9: Case 1.3 

This time we have to prove 

Mt(a, f))f (f) - Mt(a n ff +, f n ff +))f (f n ff +) 
- At(t(a n if", f n #-))i(f n if - ) + ^(t(a n ff, f))i(f n if) 

In this case, the intersection of the transverse cones is 

t(o n ff + , f) n t(o n H~, f) = t(o n ff, f), 

the union of the transverse cones is 

t(o, f) = t(o n H + , f) u t(o n ff-, f). 
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Thus we get 

(Mt(a, f)) - ^(t(a n H + , f)) - /i(t(a n H~, f)) + /i(t(a n f)))/(f). 

From the induction hypothesis applied to the space V/ lin(L), we deduce that 
this is equal to 0. 




Figure 10: Case II 

II) Next, we consider an affine subspace L such that L n a is not a face 
of a but, say, f = L n a fl iJ + is a face of a fl H + . Then we must have L C H 
so that f is a face of the three cones a fl if, a fl H + and a fl H~, but not a 
face of a. We have to show that 

/i(t(a n H + , f n # + ))/(f n # + ) + /i(t(a n h~, f n i/ _ ))/(f n H~) 
-At(t(a n f n //))/(f n //) = o, 

with 

f n H + = f n H = f n H = f. 

In this case, the union t(a fl H + , f) U t(a fl H~, f) is the projection nn n ^(a) 
of a on V"/ lin(f) ; it is not pointed, therefore, applying again the induction 
hypothesis to the space Vyiin(L), we have 

/i(t(a n H + , f)) + /i(t(a n f)) - /i(t(a n //, f)) = /iKn( f) (a)) = 

and the result follows. □ 
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Corollary 17 Let p be a polytope in V and s G Vq. Then 



M* + <M*(M)) 



veV(p) 



is equal to 1 if the point s is integral and otherwise. 

Proof. This follows immediately from the valuation property and the 
relation (see [3]) between the characteristic functions: 



where C denotes the space of linear combinations of characteristic functions 
of cones with lines. □ 

Now we show that our functions are analytic near 0. 

Proposition 18 Let a be an affine cone in V . The function fi(a) is analytic 
near 0. 

Proof. The result is true when V = (and the explicit computation 
shows that it is true also when dim^/ = 1). We will prove it by induction 
on the dimension of V. Using the valuation property, it is enough to prove 
the analyticity when a is a solid simplicial unimodular affine cone in V. Let 
v\, . . . ,Vk be primitive integral generators of the edges of dir a. If $ is a 
meromorphic function on V* such that the product 



is analytic, we denote by Res^^) the residue of $ along v\ = 0, that is to 
say the restriction to v ^ C V* of 



From the properties of the functions S and I ( Lemma ITT]) and the induction 
hypothesis, it follows that the product 



X(dir(t(p,v))) =X({0}) mod/: 



ueV(t>) 
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is analytic. Thus we want to show that Res„ 1 (/i(a)) = 0. Starting from the 
defining formula (JUJ), we want to prove that, for £ G v j 1 , we have 

(17) Res Vl (S(a))(0= E M(t(a, f))(£) R eSi;i (/(f))(£). 

f;dim f>0 

Let us denote by tt the projection V — >• VJ < V\ >. The cone 7r(a) is also a 
simplicial unimodular cone with primitive integral generators 7t(v2), ■ ■ ■ , ^{vd}- 
If a is unimodular, then the parallelepiped s + □ (t>i, . . . , f<z) C V contains 
exactly one integral point. Therefore, the explicit computations (jUJ) and (jHJ) 
of /(a) and 5(a) imply immediately that the residues along v\ — of the 
functions S'(o) and /(a) are given by: 

Res Wl (5 , (a)) = -5 , (7r(a)), 

Res wi (/(a)) = -/(tt(o)). 

In the sum (117)1 . only the faces f for which t>i is an edge of dir f contribute, and 
these faces are in one to one correspondence with the faces of ir(a). For such 
a face f, the transverse cone of 7r(o) along 7r(f) coincides with the transverse 
cone t(a, f) and we have also 

Res, 1 (/(f)) = -/(7r(f)), 

whence ()17|). and the proposition. □ 

Next we will show that Formula Q still holds when a is replaced by any 
polyhedron p. This will be an easy consequence of Brion's theorem and the 
valuation property of / and S. In the following three theorems, we collect 
the results of this section. 

Theorem 19 Let V be a rational space and Q a rational scalar product on 
V* . IfW = V/L is a rational quotient space ofV, we denote byC(W) the set 
of affine cones in W . For a G C(W), let 1(a) and S(a) be the meromorphic 
functions with rational coefficients on W* defined in Propositions^ and^ 

There exists a unique family of maps fiw, indexed by the rational quotient 
spaces W of V, such that the family enjoys the following properties: 

(a) fi w maps C(W) to H(W*), the space of analytic functions on W* , 
with rational Taylor coefficients. 

(b) IfW = {0}, then M{0}) = 1. 
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(c) If the affine cone a G C(W) contains a straight line, then Hw{o) = 0- 

(d) For any affine cone a in W , one has 

S(a)= MTy/Mf)(t(a,f))/(f) 
where the sum is over all faces of the cone a. 

As in all this section, in Formula (d), the function fi W / hn ^(t(a, f)) is consid- 
ered as a function on W* itself by means of the orthogonal projection W* — > 
(W/lin(f))* = (lin(f))- 1 with respect to the scalar product on W* C V*. 

Theorem 20 The analytic functions defined in Theorem[TQ have the follow- 
ing properties: 

(a) For any x G A, one has fi{x + a) = /u(a). 

(b) The map a i— > //(a) is equivariant with respect to lattice-preserving 
isometries. In other words, let g be an isometry of W which preserves the 
lattice A. Then ^(^(a))^" 1 ^) = //(a)(0- 

(c) If W is an orthogonal sum W = W\ © W2 and a« is an affine cone in 
Wi for i = 1, 2, then 

fi(ai + 02) = /u(ai)/i(a 2 ). 

(d) For a fixed s G Wq, the map c — > fi(s + c) is a valuation on the set of 
cones in W . 

(e) Let p C W be a polyhedron, then 

(is) s(p)(o= £ Mt(M))(0/(f)(0- 

Proof. In Theorem QUI only point (e) has not yet been proven. If v is a 
vertex of p, let us denote by jF(p,t>) the set of faces of p which contain v. 
For such a face f G jF(p,t>), the intersection c) =< f > nt(p,v) is a face of 
the cone t(p,v) and this correspondence is a bijection between jF(p,t;) and 
jF(t(p,t>)) with inverse given by f = fl p. Moreover the transverse cone 
t(p,v) of p along its face < f > flt(p,t>) coincides with the transverse cone 
t(p, f) of p along f. Therefore we have 

S(n(p,v))= J2 M(t(p,f))/(<f>nt(p,^)). 
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Replacing S(t(p,v)) with the right-hand side of this equality in Brion's for- 
mula, we obtain 

S ^= E E Mt(P,f))/(<f>nt(p^)). 

Then we reverse the order of summation and get 

S(p)= Mt(M)) E n<f>nt(p^)). 

fe^(p) oev(0 
The last sum is equal to /(f). □ 

Theorem 21 Assume V = M. d with a fixed dimension d. Then, for m fixed, 
there exists a polynomial time algorithm which computes //(a) at order m for 
any affine cone a C V. 

Proof. By [2], there exist polynomial time algorithms which compute the 
functions /(c) and 5(c) at order m for any affine cone c in M fc , if k < d. 
Therefore by induction we get an algorithm which computes /i(a) for any 
ocV.D 

Let a be a cone in the dual space V*. The dual cone a* C V contains 
the vector subspace < a > _L . Let us denote by n <a> ± the projection V — > 
V/ < a > ± . For any s G Vq, the projected cone vr <(7> x(s + cr*) is a pointed 
cone in V/ < a > ± . Thus n(n <(T> ±(s + cr*)) is an analytic function on 
(V/ < a > L )* =< a >C V*. We consider it as a function on V* by means 
of orthogonal projection, as before. We obtain a map /i* : C(V*) — > Ti(V*) 
defined by: 

Definition 22 

= Kn <(T> ±(s + o-*)). 

From the valuation behavior of /i, it follows that fi* is a solid valuation. In 
other words, the following corollary holds. 

Corollary 23 Let a be a cone in V* , and let {a'} be a subdivision of a. For 
any s G Vq, we have 

dim cr'=dim it 
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Proof. As {a'} is a subdivision of a, we have 

dim cr'=dim a dimcr'<dimcr 

Let L =< a >^C V and let ir denote the projection V — > V/L. The map 
r l— >■ x( 7r ( s + r *)) is a valuation on the set of cones in V* (see [3] for instance). 
Therefore 

x(7r ( s + (T *))= £ x(7r(s + ,/*)) + £ ± X (7T(s + 0)- 
dimcr'=dim(T dimcr'<dim(T 

If dimcr' < dim a, then the cone ir(s + a'*) C V/L contains a straight line , 
thus fi(ir(s + &'*)) = 0, and the corollary follows from the valuation property 
of /i. □ 

In a companion paper |J, we will prove the following theorem, which 
extends to equivariant homology a result of Pommersheim-Thomas [20] by 
which they answered a question of Danilov |14j . 

Theorem 24 Let S be a fan in V* and let X be the corresponding toric 
variety. For a G S, let X(o~) C X be the corresponding orbit closure. Then 
the equivariant Todd class of X is equal to 

£ /iSto [*(*)] 

in the equivariant homology ring of X . 

5 Local Euler-Maclaurin formula 

As in the previous section, V is a rational space and we fix a scalar product 
on V. Let p be a (convex rational) polyhedron in V. To each face f of p, we 
are going to associate a linear differential operator D(p, f) on V. 

To any analytic function $(£) on V*, defined near 0, there corresponds 
a unique linear differential operator -D($) (of infinite degree) with constant 
coefficients on V such that $(£) is the symbol of -D($). More precisely, for 
^ E V*, let us denote by the function x i— > e^' 2 ^ on V, then -D(<3>) is defined 
by the relation 

/}($) ■ = $(£)e^ f° r £ small enough. 
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Let W = V/L be a quotient space of V and let a be a pointed affine cone 
in W . In the previous section, we constructed an analytic function /i(a) on 
W* = L C V*. By orthogonal projection, we consider fj,(a) as a function 
on V* and we introduce the corresponding differential operator D(fi(a)) on 
V: 

(19) J D(/i(a))-e« = /i(a)(Oe 5 . 
Let p be a polyhedron in V. 

Definition 25 Let f be a face of p. We denote by 

L>(p,f) = L>(//(t(p,f))) 

the differential operator on V associated to the transverse cone t(p, f) of p 
along f. We denote its constant term by u(p,f). Thus 

i/(p,f) = J D(p,f)-l = A*(t(p,f))(0). 

The operator D(p, f), as well as its constant term z/(p, f), are local in the 
sense that they depend only on the class of t(p, f) modulo integral transla- 
tions. In particular, if p has integral vertices, then D(p,f) depends only on 
the cone of transverse feasible directions at a generic point of f . The operator 
D(p,f) involves only derivatives in directions orthogonal to the face f. 

We are now ready to state the local Euler-Maclaurin formula for any 
polytope. 

Theorem 26 (Local Euler-Maclaurin formula) 

Let p be a polytope in V. For any polynomial function h(x) on V , we 
have 

(20) E l D M- h 

where the integral on the face f is taken with respect to the Lebesgue measure 
on < f > defined by the lattice A fl lin(f) . 

Proof. 

The method is to check equality ()20|) for a polynomial of the form h(x) = 
Taking Taylor series, we may replace h(x) by e*^' 2 ^ with t small. 
Then the equality (j2*Uj) becomes the formula in Theorem EUJ (c) 

s(p)(0= E Mt(P,f))(0/ f (0- 
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□ 

In dimension 1, when p is an interval [a%, a 2 ], applying Formulas (fT3j) and 
([HI) for ii(ax + R + ) and /x(a 2 - we obtain 

— 6((n+l),ti)„d, n 



/12 
Jl(f)dt - £ 
, „., >■- -I i, -(I 



(n + 1)! 



n>0 



where t\ and t 2 in [0, 1[ are defined by t% — k± — a±, ti = a 2 — &2, with k\ the 
smallest integer greater or equal than a\ and fc 2 the largest integer smaller or 
equal than a 2 (Figure EJ). Of course, when ai and a 2 are integers, we recover 
the historical Euler-Maclaurin formula. 



ai a 2 

Figure 11: Euler-Maclaurin for the interval [ai,a 2 ] 



6 Ehrhart polynomial 

Let p be a rational polytope in a <i-dimensional rational space V and let q 
be an integer such that gp has integral vertices. Let h(x) be a polynomial 
function of degree r on V. For any integer t, we consider the dilated polytope 
tp and the corresponding sum 

s{tp,h)= h ( x )- 

xetpnA 

As a function of t, it is given by a quasipolynomial: there exist functions 
t i — ► ^(p, /i, t) on Z which are periodic with period g such that 

d+r 

(21) 5(tp,/i) = ^^(p,/i,t)f 

8=0 

whenever t is a positive integer and even in a slightly larger range including 
negative values. 
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Definition 27 The periodic functions Ei(p,h,t) defined by Equation \21)) 
are called the Ehrhart coefficients for the polytope p and the polynomial h. 

When h is the constant polynomial h(x) = 1, we denote the Ehrhart 
coefficients simply by Ei(p,t). 

Our local Euler-Maclaurin formula gives an expression of the coefficients 
Ei(p, h, t) in terms of the functions /u((t(p, f)), as we will now explain. 

Let a be a pointed affine cone with vertex s in V. If t is a nonzero 
integer, we define //(a, t) = fi(ta). For t = 0, we define /i(a, 0) = yu(dir(a)). 
Then /i(a, t + q) = /i(a, t) for any integer q G N such that the point qs is 
integral. Let p be a polyhedron and f a face of codimension m. We define 
D(p,f,t) = D(^(t(p,f),t)) for any integer t G N. 

Remark 28 // the affine span < f > of the face f contains an integral point, 
then the vertex of the transverse cone t(p, f) is integral, therefore /i(t(p, f), t) = 
/x(dir(t(p, f))) and D(p,f,t) do not depend on t . We have D(p,f,t) = D(p,f) . 

Let Wi, . . . ,w m be an integral basis of the subspace lin(f)- L<3 c V. The 
operator D(p, f,t) has the following expression: 



where A = (a x , . . . ,a m ), with a { G N and d A = D{wi) ai ■ ■■D(w m ) am . The 
coefficients z^^(p, f, t) G Q are periodic with respect to t, with period equal 
to the smallest integer such that qf < f > contains an integral point. 

Proposition 29 Let p be a rational polytope and h a polynomial function of 
degree r on V . Then, for any integer t > 0, we have 



D(p,f,t) = u (p,f,t)+ "A(P,f,t)d A 



A,\A\=1 



(22) 




Furthermore we have 




dim f+r 



where the coefficients E^p, h, f, t) are periodic with period qf. 
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Hence the Ehrhart coefficients are given by 

E l (p,h,t) = E { (P,h,f,t). 

f,dimf<i 

Proof. For t > 0, Formula (|22j) is just the Euler-Maclaurin formula of 
Theorem EH and the second equation follows from obvious estimates on the 
polynomial behaviour of the integrals. 

For t = 0, as both sides of (j2H> are quasipolynomials, they take the same 
value. We may also deduce the equality for t = from Corollary [T7| Indeed, 
for t = 0, the left hand side is h(0) and the faces of dimension > give a 
zero contribution to the right hand side. The equality becomes 

h(0)= £ (D(dir(t(M))./0(0), 

i>ev(p) 

and follows immediately from Corollary El D 

For instance, let /ibea monomial = x m with m G N d . The coefficient 
of the highest degree term t d +H is the integral f x m dx. As the operator 
D(p,p,£) is equal to 1, the face p does not contribute to the coefficients of 
lower degree. The coefficient of £ d+ M -1 involves only the faces of codimension 
1, the coefficient of t d +\ m \- 2 involves only the faces of codimension 1 and 2, 
etc. 

When we apply the last proposition to the function h(x) = 1, we obtain 
a formula for the number of integral points in tp. 

Corollary 30 (a) The Ehrhart quasipolynomial of the polytope p is given by 
Card(tpnA)= ^(p,M)vol(f)t dimf , 

hence 

E k (p,t)= ^o(p,f,t)vol(f). 

f,dim f=fc 

(b) The rational number u (p,f,t) depends only on the class modulo lattice 
translations of the transverse cone tt(p, f) . Therefore, it is a periodic function 
of t with period at most equal to qf, the smallest integer such that gj < f > 
contains integral points. In particular, if the afftne span < f > contains inte- 
gral points for every k-dimensional face f of p , then the Ehrhart coefficient 
Ek(p,t) does not depend on t. 



30 



(c) When dimV^ is fixed, there exists a polynomial time algorithm which 
computes z/ (p, f, t). 

Proof. The last statement in (b) is due to Stanley [221; it also follows imme- 
diately from Remark EHl The computability follows from the computability 
of the functions /i(a). □ 

Barvinok [3] proved recently that, given an integer m, there exists a 
polynomial time algorithm which computes the m highest coefficients of the 
Ehrhart quasipolynomial of any rational simplex in M. d , when the dimension 
d is considered as an input. We hope that our construction of the functions 
fj, leads to another polynomial time algorithm which would compute the m 
highest coefficients of the Ehrhart quasipolynomial for any simplex in M. d and 
any polynomial h(x), when the dimension d and the degree of the polynomial 
are considered as input. 



7 Computations in dimension 2 

In this section, V = R 2 with A = Z 2 . Let t>i and v 2 be primitive integral 
vectors. Let s G Vq. We are going to compute /x(a) for the affine cone 

a = s + R+^i + R+v 2 - 
We will use the following notations: 

For t G R, we denote the smallest integer greater or equal than t by ceil(t) 
and we define [[*]] G [0, 1[ by [[*]] = ceil(t) - 1. 

We denote by B the function B(y,t) = + ^- Recall that the function 
fi(a) for a one- dimensional cone (half-line) a is expressed in terms of B. 
s = SiVi + s 2 v 2 with Si G Q. 
d = g^l, for % = 1, 2. 

q = det(t>i,f 2 ). We assume q > 0. 

w G Z 2 is a vector such that det(vi,w) = 1, 

p = det(t>2, w). Thus p and g are coprime integers. 

t = (q s i + [ND+p(?sj + [WD- 

^ is a primitive g-th root of 1. 

We observe that the lattice in V/Mvi is generated by w = where w, v 2 
denotes the image of w, v 2 in VfRv%. Thus the transverse cone t(o, fi) C 
V/Mvi is given by 

t(0,fi) = (qs 2 + R + )-v 2 . 

Q 
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Proposition 31 For £ e V* , let yi = (£, v^), fori = 1,2. T/ien the function 
/i(a)(£) zs given by 

(23) /.(o)(£) = 



q-l tkr 

E 



9 \ v (i_ e T)(i- e T) (1 - • I ) 

. 1 T y ( y2-C 1 y 1 1 y x -C 2 y 2 , r| , 

+— £( — - — , [MD + — b{ — - — , [[?si]]) 



i/i <? z/2 g 2/12/2 

Its value at £ = is equal to 

(24) //(a)(0) = i((i-[[g S i]])(i-[M)) 

+ ^ " 2 [bS2]] + 2 [[qS2]] } + QR^) ( 12 - 2 [[qSl]] + 2 [M] } 

+ X ^ 



1 

a 



g^(l-C fc )(l-C fcp )' 



When a is unimodular, /i(a)(£) is given by 
exp([[si]]j/i + [[s 2 ]]y 2 ) 



(25) /i(a)(0 



(1 - e»)(l - e^ 2 ) 
+ -5(^/2 - C m , [[s 2 ]]) + -5(yi - C 2 y 2 , [[ 8l ]]) 



2/i 2/2 2/iJ/2 

an<i its constant term /i(a)(0) is given by 

(26) ^(a)(0) = (i-[[ ai ]])(i- W) 

+ QK^) ( 12 - 2 M + 2 M } + QK^) ( 12 - 2 M + 2 [[Sl]] } ' 

Remark 32 In actual Maple computations, we use only the unimodular case 
i25\) which is computable in polynomial time at any given order. Thanks to 
the valuation property, for a non unimodular cone a, we compute /i(a) by 
performing first a signed decomposition of a into unimodular cones, similar to 
Barvinok's decomposition. As a result, by our local Euler-Maclaurin formula, 
we have fast algorithms which compute, for a polygon pel 2 and a monomial 
h(x) = x^x™' 2 , the sum of values at integral points S(p,h) = J2xepnA^( x ) 
and the coefficients of the corresponding Ehrhart quasipolynomial. 
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Proof. We use the defining relation of Proposition rjjZJ First, we obtain a 
summation formula for S(a) by using finite Fourier transform as in jU]. We 
observe that AcM = ^(Zt>i + Zi^). Let a C a be the cone 

a = s + R + V! + R + v 2 

with vertex ^ 

s = -(ceil(gsi)fi + ceilfg^)^)- 

q 

Then a H M = a n M. As A C M, we have also a n A = a n A. Consider the 
dual lattice M* C A*. Let x G M. We have 



E < 

7GA* /M* 



2^(7^) = if x i Z 2 
g if x G Z 2 . 



Therefore we have 



= ; E E 

Since 



9 

7GA*/M* zG&nM 



e (2i7T7+^,x) 



an m = s + z+— + z + — 



we obtain 



1 9-1 (2iirkS,s) 

(27) g(a)(0 = g(a)(0 = e&>- £ (^^h ' 

q k=o( 1 - e c '« ; )(l-e^ c '«') 

In this formula, 5 is a generator of the group A*/M* = Z/gZ. By using the 
basis (vi,w) of Z 2 , we obtain: 

(6,vi) = l, (5,v 2 )=p. 

Let (v*,^) b e the dual basis of (fi,^)- The orthogonal projection of £ on 
(y/Kvi)* = is ec l ual to (-Ci^i + V2)v2, with 

= QK,^a) = Q(vi,V2) 
1 Q(«a,«a) Q^i^i)' 
Then the computation in dimension one (jll)) gives 
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and similarly 



with 



We have 



Co 



Q{Vl,V2) 

Q{v 2 ,v 2 ) 



/(fi)(0 = c B1 ' l4wa (— ). 

J/i 

J(a)(0 = e ^ 1+y2S2 — . 
Therefore, by (|I2J). we have 
(28) //(a)(0 = e- felSl+y2S2) 5(a)(0 

+ — , gs 2 ) + — S( , [[qsi }}) 

Vl Q 2/2 q 

q 

2/12/2 ' 

In Formula (J2EJ), we replace 5(a) with the right hand side of ([27]). taking in 
account the equality ceil(gSj) — gsj = This gives (J23j) . □ 

If a is not unimodular, then /x(a)(£) involves the "extended" Fourier- 
Dedekind sum 



^(l-C fe eV)(l-C^eV) 
and //(a)(0) involves the Fourier-Dedekind sum 

qfr! (i-C fe )(i-C fcp ) 

One has (see for instance [Sj) 

Bfe.i, P , r ) = £((-^))(A)-i-, 

fc=0 H H H 

where the '"sawtooth" function ((a)) is defined by 

((a)) = a - [a] - ~. 
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The Dedekind sum D(q, l,p,r) can also be computed in polynomial time by 
means of reciprocity relations (see for example but here we do not use 
this fact. 

5131761430387/12155092992 

287696501/133706022912 

1/10626 

28224572717107/66853011456 



-1/252 

Figure 12: Le savant Cosinus 



Example 33 ( Fiaur iT^) We compute the right hand side of Euler-Maclaurin 
formula in the case of the "dull triangle" with vertices (0, 0),(1, 0),(0, 1), and 
the polynomial h(x) = xfx 2 - As expected, the contributions of the various 
faces of p add up to 1 . 

Contribution of vertices- (4 28224572717107 5131761430387 
^Ulbbl lUULtvib UJ ueiiices. V, 66853011456 ' 12155092992 ' 

Contribution of edges: — ^y i^^qI^u , 0. 
Integral over triangle: 1Q q 26 . 



Example 34 Triangle with vertices Si = (|, |) ; s 2 = (y, \), s 3 = (y, y). 
Number of integral points: 31 

Contribution of vertices- 89133678169939 4281800310619 401172431621091 
<^unil luuiivn uj UKibiLKd. 660 88208614500 ' 2106396270216 ' 457987274773000 

Contribution of edges /yy — yy y^Q 
Area of triangle: y|y 

Example 35 Quadrangle with vertices (~, |) ; (y, |), (y, y), (3, 10). 
Number of integral points: 49 

Contribution of vertices- 210849514883 4281800310619 179008247 4382929 
kjuiiui uuu,ui,uii uj uci uul,co. 127 95 6 3 2 2980 ' 2106396270216' 706816180' 6869864 

Contribution of edges : 2T/j ; ~ yo' M' M 
Area: y^ 

1 This computation delighted us, and it would have delighted Dr. Pancrace Eusebe 
Zephyrin Brioche alias "Dr. Cosinus" ^Sj 
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401172431621091/ 457987274773000 



89133678169939/6608820861450C 



-4382929/68€9SM 



210849514883/12795632298C 




T/210 



4-2&1800310619/2106396270216 



179q08247/ i706816180 



T/210 



4:2&1800310619/2106396270216 



Figure 13: Triangle with vertices (§, |), (f , |), (f , 92 
extra vertex (3,10) 



^) and quadrangle with 



Remark that, as expected, the contributions of the bottom right vertex (y , |) 
in the triangle or the trapezoid of Figure I*HI1 are the same, as this vertex have 
the same tangent cone in both polygons. 

Example 36 We compute the Ehrhart quasipolynomial E 2 t 2 + Ei(t)t + E (t) 
for the number of integral points of the triangle of Example \34\ The highest 
coefficient is the area of the triangle, E 2 = . The coefficient E 1 (t) is the 
sum of the contributions of the edges. The coefficient E (t) is the sum of the 
contributions of the vertices. 

On this example, we can observe the periods of the contributions of the 
edges and vertices to the Ehrhart coefficients ( Corollaru \3U\) . The period of a 
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vertex contribution is equal to the 1cm of the denominators of its coordinates. 
For an edge starting from a vertex (a±, a 2 ) and parallel to the primitive vector 
(vi,V2), the period is the least integer q such that q(a\V2 — a 2 vi) is an integer. 

Contribution of edges (periods 3, 5, 7 respectively) : 

— 155 mod (f,3) + jq, mod (4 1, 7) + 1/30, mod (2 1, 5) + ^ . 

Contribution of vertex si = (|, |) (period 15): 

^ mod (t, 5) mod (t, 15) — ^ mod (t, 3) mod (t, 15) — ^ mod (t, 5) mod (2 1, 3) 

+ 1 BS - T ™ rf (*» 5 ) - ilia ^ ( 2 *. 5) + I ^ (7*, 15) 
+ | mod (8 t, 15) - i mod (i, 15) + § (mod (t, 5)) 2 

(mod (8 t, 15)) 2 + (mod (i, 15)) 2 + § mod (2 i, 3) 

+ ^mod(t,3)-| (/mod (2 t,3)) 2 + ff| (mod(t,3)) 2 

" ngsSo (mod(2t,5)) 2 - ife (™<>d (7t, 15)) 2 + ^ mod (2 1, 5) mod (7t,15) 
+ ^ mod (8 i, 15) mod (2 1, 3) 



Contribution of vertex s 2 = (y, |) (period 21): 

- 477S776 5 - I ( morf (*> 7 )) 2 + f ™ orf (*> 7 ) + es" ™ orf (*> 3 ) ( 2 *> 21) 

mod (t, 3) - i mod (2 i, 21) - ^ (mod (t, 3)) 2 + § (mod (2 i, 3)) 2 

-§ mod (2 i, 3) + i mod (2 i, 21) mod (2 i, 3) 

+ ^mod (2t,3) mod (4t,7) - mod (At, 7) 

-SM(^.7)) 2 
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Contribution of vertex s 3 = (y, y) (period 35): 

(mod(W t, 35)) 2 + ^ (mod(23 f , 35)) 2 + gSSo 
+ JH mod(4*, 7) + ^ (mod(3At, 35)) 2 
+l/10mod(3t,5) - 1/35 mod(34t,35) - ^ (mod(18 1, 35)) 2 
+ _|_ ( mod (g f ? 35 ))2 + ( mod{2 1, 5 )) 2 

+ ^ mod(18 *, 35) + ^ (mod(3 t, 5)) 2 + ^ mod(34 1, 35)mod(18 1, 35) 
+ j^fmod(23t, 35)mod(16t, 35) 

+ y|^ mod(3t, 5)mod(23t, 35) — mod(34t, 35)moe*(16t, 35) 

— -p^ mod(2t, 5)mod(16t, 35) 

— 1/35 mod{3 1, 7)mod(3 1, 5) 

— 1/35 mod(3t, 5)mod(At, 7) — mod(9t, 35)mod(3t, 7) 

— ^5 mod(3t, 5)mod(9t, 35) 

mod(16 1, 35)mod(18 t, 35) + yggff^ mod (2 t, 5) 
+ s^fes ( moc? (4 7 )) 2 - 1/35 mod(9 1, 35) - ^ mod(23 1, 35) 
-J mod(16 1, 35) + 1/49 (mod(3 1, 7)) 2 

Example 37 (Computation time). 

We computed the full Ehrhart quasipolynomial corresponding to the tri- 
angle of Example |3"""| and the polynomial h(xi,x 2 ) = x^x^ 2 , with increasing 
exponents k\ and k 2 . Allowing a computation time of about one hour, we 
reached k\ = k 2 = 24. The result is of course too big to write here. 

The sum of values xf'x^ at the integral points of the triangle of Example 
pm dilated by the factor N = ll 5 took about the same time. The result is the 
following number 

55969247458735493271268368615238071121335974262337882261418363621 
89704055956429496253759473056373507451253522021344188115187647607 
84555431172202923756940824265247663088847763429436570335188702325 
06644969965841257822711805056447218921550669146263582661876630783 
21357671611262065293901983868557252464459832189159990869820527095 
53646871654914800005753059422066576204781923454823934475242960034 
42199041253798398004263030681714027295470241663946228744550160085 
43856624239377702107746492579014275563017167813144052693763385569 
75239252588060279466314599314734680953729093269435217987689840619 
0740089242444014302. 

As experiments showed, our method for the computation of J2xepnz d h( x ) 
is very efficient for this small dimension, compared to other available soft- 
wares. Furthermore, as Example 13*61 shows, the Ehrhart polynomial is written 



38 



as a sum of canonical contributions of all faces, once the scalar product is 
fixed. We will come back soon to the computational and complexity aspects 
of this problem for higher dimensions. 
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